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In this paper the idea of holonomic quantum computation is realized within 
£N) \ quantum optics. In a non- linear Kerr medium the degenerate states of laser 

beams are interpreted as qubits. Displacing devices, squeezing devices and 
' interferometers provide the classical control parameter space where the adia- 

batic loops are performed. This results into logical gates acting on the states 
q ■ of the combined degenerate subspaces of the lasers, producing any one qubit 

CN , rotations and interactions between any two qubits. Issues such as universality, 

\ complexity and scalability are addressed and several steps are taken towards 

Q\ ' the physical implementation of this model. 

■ I. INTRODUCTION 

a ■ 
a ■ 

Holonomic transformations have been recently proposed |l|] and more extensively studied 
H as logical gates for quantum computation The idea formal as it may apear at a first 
glance is not confined to a purely theoretical shere, but has the challenging possibility of 
^ . experimental implementation. Towards this purpose we employ existing devices of quantum 
optics, such as displacing and squeezing devices and interferometers acting on laser beams 
in a non-linear medium. A different setting of optical quantum computer has been reported 
in ||]. The attempt to apply the abstract idea of holonomic quantum computation (HQC) 
to a physical system has lead us to deal with and clarify some theoretical problems of HQC 
such as universality, complexity (tensor product structure of qubits) and scalability. On the 
other hand, the experimental setup of the scheme proposed may prove to be a possible even 
though challenging task for the experimenters. 

The basic idea of HQC is related with the geometrical phases generated by the isospec- 
tral transformations of an n-fold degenerate Hamiltonian, Hq, so its presentation can be given 
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in a geometrical form |T| . Initially, quantum information is encoded in the n dimensional de- 
generate eigenspace C of H , with eigenvalue E . The operator H Q is considered to belong to 
the family T = {H a = U(a) H W(a); a G M} of Hamiltonians unitarily (U^a) = W l {a)) 
equivalent and therefore isospectral with H , where H = H ao for some a G M.. as a ranges 
over the control manifold Ai no energy level crossing occurs. The <r's represent the classical 
"control" parameters that one uses in order to manipulate the encoded states G C. Let 
C be a loop in the control manifold M.. When the loop C is slowly gonne through, then 
no transition among different energy levels occurs and the evolution is adiabatic, i.e. T is 
faithfully realized by the experimental setup. If \ip)i n G C is an initial state in the degenerate 
eigenspace, at the end of the loop it becomes \ip) ou t = e tE ° T Ta(C) \ip)i n . The first factor is 
just an overall dynamical phase which in the following will be omitted by a redefinition of the 
energy levels, taking Eq = 0. The second contribution is the holonomy Ta{C) G U(n), and 
is a result of the non-trivial topology of the bundle of eigenspaces over M.. By introducing 
the Wilczek-Zee connection [[J 

AZ := (p\U\a) ^-U{a)\p) , (1) 

where A^f is the (p, p) matrix element of the <7j component of the connection, one finds 
T A (C) = Pexp/ C A, [§, where P denotes path ordering. The set Eo\(A) := {T A (C);VC G 
Ai} C U(n) is known as the holonomy group 0. In the case where it coincides with the 
whole unitary group U(n) the connection A is called irreducible The transformations 
Ta(C) for suitable C's can be used as logical gates for the HQC. 

We shall focus on quantum optics, a well established area of quantum physics, in which 
the developed technology is quite mature as a possible venue for practical implementation of 
HQC. The model we study here includes laser beams moving through non-linear Kerr media, 
and acted on by displacing and squeezing devices and interferometers. This implementation 
has the merit that it gives direct answers to several problems which were raised in the 
theoretical study of HQC 0. 

In Chapter II we present the schematic theoretical description of the quantum optical 
components employed for HQC. This includes the non-linear Kerr medium, the one and two 
mode displacing and squeezing devices as well as their effect on the states of laser beams. 
In Chapter III we construct the non-Abelian Berry connection, the field strength and the 
holonomies related with this optical setup. A model with SU(2) interferometers is also 
given as an alternative tool for classical control, and its holonomies are calculated resorting 
to the non-Abelian Stokes theorem. In Chapter IV the connection between the experimental 
components of quantum optics and the theoretical requirements for HQC is described. A 
numerical simulation is finally reported indicating the reliability of the logical gates with 
respect to the scale resources of the HQC. In the Conclusions the quantum computation 
characteristics of this model are discussed and issues like the universality, complexity and 
scalability are addressed. 

II. THE QUANTUM OPTICAL MODEL 

In the following we shall exploit the advanced tools of quantum optics in order to imple- 
ment a specific HQC model. All the components used here are thoroughly analyzed in the 
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optics literature || and experimentally realized by employing such devices as beam splitters, 
frequency converters, four wave mixers, and others. 

A. Kerr medium Hamiltonian and Degenerate States 

In order to perform holonomic computation we shall employ the nonlinear interaction 
Hamiltonian produced by a Kerr medium 

Hj = hXn(n - 1) , 

with n = a^a the number operator, a and being the usual bosonic annihilation and 
creation operators respectively, and X a constant proportional to the third order nonlinear 
susceptibility, x ■> °f t ne medium. Degenerate eigenstates of Hi are |0) and |1) ({|z/);z/ = 
0, 1, ...} denoting the Fock basis of number eigenstates n\v) = v^))- In the case of two laser 
beams, with annihilation operators ai and a 2 respectively, the total Hamiltonian is given by 
the sum 

H} 2 = hXmint - 1) + hXn 2 (n 2 - 1) . 

Its degenerate eigenstates are the tensor product of the eigenstates of each subsystem: \i1j2) 
:= (g> lj'2) for H,j 2 = 0, 1 with |ix) and \j 2 ) the degenerate states of each beam. Accord- 
ingly, the unitary transformations acting on the system are given by the tensor product of 
the transformations on each individual subsystem. For example, the transformation of a 
system (Hamiltonian and states) of two lasers when one beam is transformed by U\ is given 
by the tensor product Uu — U\ (g) 1. These rules can be applied to build up a system with 
m lasers. In this case the subspace of Fock states on which we restrict in order to apply the 
adiabaticity theorem has as basis vectors the degenerate states |0/) and for each laser la- 
belled by I. The general state of the system of m lasers is given by \p\...p m ) = \p\) \p m ) 
where pi could be zero or one, for I = 1, m. On this space of states the code can be writ- 
ten. We have good reasons to believe that the problem of the generation of stable Fock 
states will be overcome, as suggested by some recent developments Q. 

B. One and Two Laser-Qubit Transformations 

On state of a laser beam with annihilation operator a, the following operators can 

act 

Displacer: = exp(Aa^ — Aa) , 

where A is an arbitrary complex parameter. The displacing device that implements -D(A) is 
a simple device that performs a linear amplification to the light field components. 

Squeezer: S(p) = exp(/ia^ 2 — pa 2 ) , 

where p is an arbitrary complex parameter. The squeezing operator S(p) can be imple- 
mented in the laboratory by a degenerate parametric amplifier. 
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The transformation operators D(X) and S(fi) acting on a single laser beam will result, 
after a closed loop is performed in their parameter space, into rotations in the state space 
spanned by |0) and |1), according to the adiabatic theorem. 

The displacer D(X), transforms the operators a, and any analytic function thereof 
f(a, a'), for any choice of parameters A, as follows QTUJ 

D(X)aD^X) = a - X , D(X)aW^(X) = o) - X , 



D{X)f{a, at) J D t (A) = f(a - A, a 1 " - A) . 
Similarly for the squeezing operator 

S(fi)aS^(fi) = cosh(2r)a + e~ %e sinh(2r)a^ 



S(n)dtS*(n) = e ie sinh(2r)a + cosh(2r)at , 
S(fx)f(a, aJ)S^(fj,) = f (S(/j,)aS^(iJ,), S(fx)a jl S^(fx) 



where /x = re 10 , with r > and —tt < 6 < tt. 

On the general state of two lasers 1^12) = l^i) <8> \ip2) with corresponding annihilation 
operators a\ and a-i-, the following operators can act 

Two mode squeezer: M(() = exp((a{a2 — C a i a 2) ■ 

The operator M((), can be implemented in the laboratory by a non-degenerate parametric 
amplifier. 

Two mode displacer: = exp(£a|a2 — ^aial) ■ 

M{C) and iV(£) are the transformations between two laser beams that produce, after per- 
forming adiabatically a loop in their parametric space, coherent transformations in the two 
qubit state space spanned by 1 00) , 1 01) , 1 10) and |11). 

These transformations on the states of the laser beams can be produced by SU(2) or 



SU(1, 1) interferometers ||11|| , according to the algebra which their generators belong to. For 
instance, each one of Xa^ — Xa and £a|a 2 — ^aia\ belongs to an su(2) algebra, while ^a) 2 ' — j2a 2 
and C Q i a 2 — C a i a 2 belong into (different) su(l, 1) algebras jl2| . 



III. APPLICATION TO THE HOLONOMIC THEORY 



The non-Abelian Berry connection, A, is generated by the topological structure of the 
bundle of the degenerate sub-spaces. It determines the way to perform a parallel transport 
of the degenerate eigenstates along an adiabatically spanned loop. In this section we shall 
show that a complete set of holonomies of A can be explicitly calculated for our model. 
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A. The Connection A 



We initially perform the following polar decomposition of the control variables 

A = r e ie ° , // = r ie i91 , C = r 2 e i(>2 , f = r 3 e id3 . 

We obtain the connection, A, from (fj), parametrizing the control manifold by the set of 
real variables introduced above M. := {r^,^; i = 0, ...,3} with elements Oi G Ai, where 
we take U(o) = D(X)S(/j i ) for the one laser transformations and U(a) = N(£)M(Q for 
transformations between two lasers. We have 



A 



-(e~ i6 ° cosh 2r x - e i{e ° +6l) sinh 2r x ) 

e ie ° cosh 2ri - e ^ e ° +dl) sinh 2r x 



A, 



ir 



[e~ Wo cosh2rx + e i( - 9o+dl) sinh2r^ 



ir (e ido cosh2ri + e ~ i{e ° +dl) sinh2r!) 



in 



For the connection components, A ro and Ag , it is more convenient to use for the variables 
A the decomposition A = x + iy, with x and y real, resulting into the following components 
of the connection 



4 = cosM r „-^\ 



'•() 



— iy —(cosh 2ri — e ldl sinh 2r\ / 

cosh 2ri — e~ ldl sinh 2r\ —iy 



A v = sm9 A ro + zz&A eo 
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The components Ag 2 and Ag z have more complicated forms that we shall not give explicitly 
here as they are not necessary for performing universal quantum computation |13| . 



B. The Commutators and the Field Strengths F 

In order to be able to calculate the holonomies it is convenient to consider loops 
C on the planes (<Ji,<jj) in -M,[] on which the two components of the connection commute 



Note that the components (ro,#o) have been replaced by (x,y). 



4 



with each other yet giving a non-trivial holonomy (i.e. they have non-zero field strength 
component, F u%aj = d ai A a . - d aj A ai + [A a ., A aj \). Indeed, for a u i = 1,2,3, denoting the 
Pauli matrices, we have 

[A x ,A ri } = with F xri \ 9i=0 = -2io- 2 exp(-2ri) , 
[Ay, A ri ] = with F yn \ 6i=0 = -2ia x exp(2ri) , 
[A ri , Ag^ = with F ri g 1 = — zs 3 sinh4ri , 



[A r2 , A r3 ) = with F, 



r2r 3 le 2 =e 3 =0 



-2ia^ sinh 2r 



2 • 



[A r2 , A ra ) = with F, 



T2V3 le 2 =0,6» 3 =37r/2 



—2ia[ sinh 2r 2 



where 
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The above conditions that are satisfied on the planes (x,ri) ei=0 , (y,ri) 0i=o , (ri,9i), (r 2 , r 3 ) 02=e 
and (r 2 , r 3 ) 02=o g 3=37r / 2 allow for the explicit calculation of the holonomies for paths restricted 
on such planes. 

C. The Holonomies T A (C) 

In order to perform universal quantum computation it is necessary to produce at least 
two independent unitary gates |TJ| . In the following we shall present holonomic gates, which 
involve (any) one qubit rotations and a special class of (any) two qubit transformations. In 
detail we have 



Cj E (x,ri) 0i=o gives T A (d) = exp-i^E/ with E 7 := J^^dxdr^e 2ri , 
C n e (y,ri) 8i=0 gives Y A {C H ) = exp -ia 2 E n with En := ^ Cu) dydr x 2e 2Tl 



Cm e (ri,0i) g ives T A {C HI ) = exp-ia 3 E in with £ m := Jz iCiii) dr 1 d6 1 smb. 4n 



Civ e (r 2 ,r 3 



> ' csy e 2 =e 3 =o 



gives r A (C/ V ') = exp -io-^S/y with X/y : = J Ti ^ Civ ~ ) dr 2 dr 3 2 sinh2r 2 



Cy e (r 2 ,r 3 ) d2=0 9a=37T/2 gives r A (Cy) = exp -iaf E v with £y := J s{Cv) dr 2 dr 3 2 sinh 2r 2 , 
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where E(C P ) with p = I,...,V is the surface on the relevant submanifold (crj, aj) of M. 
whose boundary is the path C p . The hyperbolic functions in these integrals stem out of 
the geometry of the su(l,l) manifold associated with the relative control submanifold. The 
r^(C)'s thus generated belong either in the U(2) or Z7 (4) group. Considering the tensor 
product structure of our system these rotations represent in the 2 m space of m qubits 
respectively single qubit rotations and two qubit interactions, thus resulting into a universal 
set of logical gates. Their explicit constructions are similar to those presented in D for the 
CP n model. 



D. The SU{2) Control Manifold 

In what follows we discuss the employment of SU(2) interferometer as control devices 
1 1] for producing holonomies. For ax and a 2 the annihilation operator of two different laser 



beams, consider the Hermitian operators 



J X 



a\a 2 



J a 



;(a|a 2 



J, 



■(a\ai 



a\a 2 ) 



(2) 



and 



iV = a\ai + a 2 a 2 



ni + n 2 



The operators (|2]) satisfy the commutation relations for the Lie algebra of SU(2)\ [J x , J y ] = 
iJzi [Jy, Jz] = iJxi [Jz, Jx] = iJy The operator N, which is proportional to the free Hamil- 
tonian of two laser beams, commutes with all of the J's. On the other hand, however, the 
Kerr Hamiltonian does not commute with the J's, allowing for the possibility that SU (2) 
interferometers be used as transformation controllers in view of the holonomic computation. 

From these operators we obtain the unitaries, U x {a) = exp(iaJ x ), U y (/3) = exp(i/3J y ) 
and U z (j) = exp(27J 2 ). For the degenerate state space of two laser beams spanned by \i1j2), 
we have from (|l|) and for IA = U x (a)U y (P)U z ( r f) the following connection components 
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(3) 



These components do not commute with each other, when projected on planes with non- 
trivial field strength. Hence, it is not possible to employ again the method used in the 
previous section to calculate the holonomies of paths in the three dimensional control pa- 
rameter space, (a,j3, 7). Instead, for this purpose we may employ the non-Abelian Stokes 
theorem 
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The extra limitation, now, for the choice of the path comes from the con- 
straint that, apart from being confined on a special two dimensional subspace, it has to have 
the shape of an orthogonal parallelogram with two sides lying along the coordinate axis. 
This will facilitate the extraction of an analytic result from the Stokes theorem. Though, 
experimentally, this restriction poses additional (possibly minor) difficulties, theoretically it 
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leads to the very interesting possibility of a direct calculation of non-Abelian holonomies 
without resorting to their Abelian substructures. Note that the application of the Stokes 
theorem for the evaluation of the holonomies in the previous section gives the same results, 
as expected. 

To state the non-Abelian Stokes theorem let us first present some preliminaries, where 
a few simplifications are introduced, as its general form will not be necessary in the present 
work. 
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FIG. 1. The loop C for the non-Abelian Stokes theorem. 



Consider the Wilson loop (holonomy), W = P exp § c A, of the loop C given in Fig. [I], 
with connection A, made out of the Wilson lines Wi for i — 1, ...,4, as W = W4W3W2W1. 
(er, r) is a reparametrization of the plane where the loop C lies. Define T _1 (cr, r) = W4W3. 
Then, for F ar the field strength of the connection A on the plane (a, r), W is given in terms 
of a surface integral 

w = p^ e / E T- 1 ( ( 7,T)F CTT ( ( 7,T)T(<7,r) ( i<7dr 

where P T is the path ordered symbol with respect only to the r variable, contrary to the 
usual path ordering symbol P, which is with respect to both variables, a and r. Here 
F aT (a, r) = —d a A T + d T A a + [A a , A T ). 

From the connection given in (|3|) the following holonomies are derived. For a closed 
rectangular loop G\ G (a, /3)-plane with coordinates (0,0), (a = 7r,0), (a = ir,/3), (0, 0) we 
obtain the following unitary transformation 

iVd) = exp(-z2/3<r 2 12 ) . (4) 

In addition for a rectangular loop C 2 E (a,7)-plane with coordinates (0,0), (a = 7r,0), (a — 
Ti, 7), (0, 7) we obtain the holonomy 

T A (C 2 ) = ex P H27<x 3 12 ) , (5) 



where the matrix a\ is defined similarly to a\ arid a\ in Subsection [III B| . These operations 



can be implemented by using interferometers between any two laser beams. Note that the 
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coefficients in front of the matrices in the unitaries are areas on spheres spanned by the 
angles a and f3 or a and 7. This is consistent with the geometry of SU(2). 

These two matrices can produce any unitary transformation of one qubit encoded in 
a sub-space of states of the two laser beams spanned by {|01), |10)}. In other words, we 
need two laser beams to encode one qubit, contrary to previous construction. As these 
transformations can be performed between any two beams, we can generate interaction 
transformations between two qubits, resulting finally (together with the one qubit rotations) 
into a universal set of transformations. For example the SWAP two qubit gate given by 



UsWAP 



10 
10 
10 
1 



is achieved as follows. On four arbitrary laser beams 1, 2, 3 and 4 with {|01), 1 10) }i,2 encoding 



the one qubit and {|01), 1 10) }3 > 4 encoding the other, we may act with Ya{C\ 



between 



beams 1 and 3 and with Ya{C\)\ 3= ^h between 2 and 4 producing eventually the Uswap 

4 

gate. The loop G\ is defined as in (f|). 

This model facilitates the physical implementation as it will be seen in the following 
section. 



IV. TOWARDS EXPERIMENTAL IMPLEMENTATION 

We address here the task of combining the theoretical requirements of HQC together 
with the features of the "experimental" components described in the previous two sections. 
While the Abelian holonomies have been produced in the laboratory by various means, the 
non-Abelian ones are more complicated. However, the holonomies calculated above, require 
successive restrictions on two dimensional planes of the control parameter space, quite in 
the same way as one needs to do to generate Abelian Berry phases. This constructive 
method may prove experimentally advantageous for performing and measuring non-Abelian 
holonomies. A survey over some Berry phase experiments in optics is given below. 



A. Various Abelian Berry Phase Setups 

Photons can be seen as massless spin-1 bosons. This characteristic has been the driv- 
ing force for the optical manifestation of the Berry phase with respect to the polarization 
quantum numbers jllj. Necessary condition for the generation of this phase factor is the 
adiabatic change of the direction of the photon propagation. Various optical experiments 
have been performed. Results at the classical level have been reported in |17|j , for the case 



of a single mode in a wounded optical fiber, whereas quantum mechanically, in ||18|| , the 
case of a single photon has been treated. Of special importance, for our case, is the latter 
experiment where the Berry phase has been observed at quantum optical level. In this case 
the incident light is prepared in an entanglent state 

A{E')\n) E ,\n) E -E>dE' , 



S 



where A(E') = A(E — E') is the complex probability amplitude for finding one photon with 
an energy E' (\n — 1) e 1 ) or with an energy E — E'(\n = l) e-e 1 )- This type of states can be 
produced in the lab by driving a single-mode ultraviolet laser into a nonlinear optical 
crystal. A Michelson interferometer has been used for the observation of the phase in the 
output state. It was found that the output state (photons in essentially n—1 Fock states) 
had an extra phase factor due to the optical-path-length difference AL of the interferometer 
plus the contribution of the Berry phase. The form of such state is given by 



I out 



-±= J A(E')\n)E'\n) E -E'{l + exp[i0(£ - E')]}dE' 



where <p(E — E') = 2irAL/ Xe-e' + (pBerry, with 4>Berry the geometrical phase predicted 
theoretically. 

Recently, an alternative approach to the geometric phase has been considered [f[9|l , through 
squeezed states of photons. Squeezed states have been found considerably interesting in the 
field of quantum optics for various reasons, as for example, the noise reduction which is 
necessary for practical applications with noise sensitivity. 

Displacement and squeezing give different contributions to the Berry phase of the Fock 
states \v). For the case of squeezing one finds that this contribution is given by 



2n+l 



Berry 



4 J 



f (cosh4ri — 1) dr\ 



Such Berry phase agrees with the form of the diagonal connection A$ 1 in Subsection [III A 



as it was to be expected. On the other hand if we perform a loop in the control parameters 
of the displacing device we expect the following Berry phase to arise 

^Berry = j (V dx ~ xd V) ■ 

The equivalent connection of displacing in the Kerr medium (A x and A y with r x = 0, in 
Subsection |111 A| ) are non-diagonal matrices, whose holonomy cannot be calculated easily. 
In fact, as it is observed by the numerical simulations in the following, the phase factors 
produced in front of |0) and |1) are not equal, due to the off-diagonal elements of A x and 
A y . This effect is related with the degeneracy structure of each model. 



B. Free Hamiltonian and Kerr Medium 

In the previous sections we used the Kerr non-linear Hamiltonian in order to produce 
the degenerate eigenspace spanned by |0) and |1). The full Hamiltonian of the system is the 
combined one of the free photons and the non-linear medium, i.e. Hxot = Hpree + Hxerr = 
hum + hXn(n — 1). Of course the first part lifts the degeneracy of |0) and |1) destroying 
the basic requirements for the holonomic computation. In order to overcome this problem 
we resort to the following constructions. 

Considering Hp r ee as unperturbed Hamiltonian and the non-linear part as the interaction 
term we may move to the interaction picture of the full system, with Hj = H^err- The 
rotation to the interaction picture may be incorporated in the devices used for the external 
control resulting in a redefinition of their control parameters. 
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Alternatively, we may define a one dimensional lattice with points on the trajectory 
of the laser. As the free Hamiltonian is acting only on the state |1) changing its phase by 
e -tH Free At _ e -ihu)Ax/c^ c the speed of light, we may single out the points Xk = 2-nck/ [Tluj) 
for k integers. On these points the phase is trivial and it does not contribute to the state. 
Hence, |0) and |1) are degenerate on this lattice |2C||. 



In Subsection |III D| we have introduced SU(2) interferometers as control devices. The 
su(2) operators commute with Hp ree allowing the effect of the free Hamiltonian to factorize 
out of the whole control procedure. At the end of the algorithm the detectors may be placed 
on a point of the degenerate lattice in order to avoid the dynamical phase produced by 
H Free on the states |0) and |1). Even though in the SU{2) model each qubit is encoded with 
the help of two laser beams increasing in this way the necessary resources, it overcomes the 
problem of the degeneracy in the most efficient way. 



C. Holonomies and Devices 

For the implementation of the continuous adiabatic loops we should adopt the kick 
method described in 0, |21] and p2| . A general state \ip) in the degenerate eigenspace of 



H = Hxerr is given as a linear combination of |0) and |1). Under an isospectral cyclic 
evolution of the Hamiltonian in the family !F, the evolution operator acting on is given 
by the 2x2 submatrix in the upper left corner of 



U(0,T) = Texp-i F U(a(t))H UHa{t))dt 
Jo 



This evolution takes place from time to time T and, for performing a closed loop, we 
demand er(0) = a(T). By dividing the time interval, [0, T], into m equal segments At we 
may approximate the above operator by 

m 

U(0, T) KTH^e-'^W] with U, = U{p>) = U{<r(U)) . 

Assuming the evolutions Uj +1 Ui to be a very small rotation and restricting to evolutions 
which remain in the zeroth degenerate eigenspace we might once more derive the holonomy 
operator Ta(C) for A defined in (0). We prefer instead to see what the effect of finitely 
many devices would be, when acting on the space of states of the qubits (the lasers). 

For the sake of concreteness we work out examples in terms of displacing devices -D(A), 
performing a closed loop in their control parameters A. This is shown in Fig. ||, where for sim- 
plicity the least possible number of displacing devices (three) for performing a closed loop has 
been considered. Two displacing unitaries are combined as _D(A)-D(A') = exp (i c 3(\\'))D(\ + 
A'). The physical process behind this is as follows. On the state first acts a displacing 
device with unitary D^(Xi), taking it to the point Ai. Then, the evolution operator of the 
Kerr Hamiltonian acts for a time interval At = T/3 U(At) = exp(—iH At). This effect is 
achieved by propagating the beam inside a Kerr medium. Then, the evolution D j; (X2)D(Xi) 
is performed. This is achieved, with a single displacing device, given (up to an overall phase 
factor that will cancel at the end) by D(Xi — X 2 ). After exiting the displacing device (we 
are at point A2) the beam enters a Kerr medium for time At and then the procedure is re- 
peated until we come back to the point Ai and the beam enters once more the Kerr medium. 
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Finally, the state is thus displaced by D(Xi). This loop may be transported to any other 
place of the control parameter complex plane by acting at the beginning and at the end of 
this procedure with the appropriate displacing unitary (device). In this case the evolution 
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FIG. 2. The triangular (and polygonal) loop C on the complex plane of the displacing control 
parameters, A, approximating the circle. 

operator is approximated by 

U(0, T) « D(Xi) (U(At; 0)U(At; A x + \ 2 )U(At; Xi)U(At; 0)) D f (Ai) , 

where U(At; A) = D{\)U{At)D^(X), \ = - \ { and A 4 = A x . 

According to the above analysis, we proceed to the numerical simulation of a system 
with various numbers of displacing devices represented by different polygons on the control 
complex plane (see Fig. |2|). We start with a pentagon which demands five displacers. Fig. 
|3] represents the absolute values of the (0,0), (0,1), (1,0), (1,1) elements of the evolution 
operator U(0, T) as functions of the number of displacers used to approximate a cyclic 
evolution. These are the relevant elements for the evolution of the states in the degenerate 
eigenspace describing a qubit. The parameters involved are taken to be T = 0.1 and hX = 1, 
with the radius of the circle equal to 1. The initial point is taken to be the origin of the 
complex plane rather than Ai, or in other words we do not perform the initial and final 
displacings by D^(Ai) and -D(Ai). In the table below are depicted the percent deviations 
of those values obtained with 5, 10, 20 and 26 displacers with respect to the ones obtained 
with 100 displacers. 





5 


10 


20 


26 


00 


0.2419 % 


0.0595 % 


0.0149 % 


0.0099 % 


01 


0.9119 % 


0.2260 % 


0.0558 % 


0.0186 % 


10 


0.9119 % 


0.2260 % 


0.0558 % 


0.0186 % 


11 


1.6763 % 


0.4061 % 


0.0760 % 


0.0269 % 



We see that with 26 displacers the error is of the order of 1 in 10 4 acceptable for quantum 
computation with error correction. This provides an indication for the necessary number of 
devices needed in order to reproduce faithfully the holonomic adiabatic loop. 
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0.398 
0.397 
0.396 
0.395 
0.394 
0.393 



10 15 20 25 30 35 



10 15 20 25 30 35 



0.398 
0.397 
0.396 
0.395 
0.394 
0.393 



0.63 



0.625 



0.62 



10 15 20 25 30 35 



0.615 



5 10 15 20 25 30 35 



FIG. 3. The absolute values |C/oo(0,T)|, \Uoi(0, T)\, |[/io(0,T)| and \Un(0, T)\ as functions of 
the number of sides of the polygons. 



V. CONCLUSIONS 

The implementation of HQC in the frame of quantum optics has provided novel insight 
into many technical aspects of the theory. Moreover, the components demanded for it are 
widely used in the laboratories. The possibility of overcoming the difficulties in combining 
them in the appropriate way for obtaining holonomies is an open problem to be faced by 
experimenters. 

In summary the main quantum computational features we observed in our model are 
the following. First, the universality condition is proven explicitly, stemming out of the 
ability to construct holonomies representing any possible logical gate. This is achieved by 
combining one qubit rotations (realized by displacing and squeezing devices) and two qubit 
transformations (by interferometers) between any two qubits. Second, the setup exhibits 
quantum entanglement, having built in tensor product structure as it consists of a multi 
partite system. This resolves the problem of complexity posed in which is one of the 
main features which make QC's more efficient than classical ones. Third, the degenerate 
space of the Hamiltonian eigenstates, which is used to write the code is constructed out of 
laser beams each with a two dimensional degenerate space. So the demand of using a big 
degenerate space to write useful codes is performed not by resorting to one system with very 
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large degeneracy, which is almost impossible to realize in nature, but by adding up the 2- 
dimensional subspaces of the lasers. This is the characteristic of scalability of the proposed 
model. Fourth, the chosen loops associated with the given holonomies are restricted on 
specific planes (<Tj,<Tj) of two control parameters Oi and <jj, exactly in the same way as 
used for the production of Abelian Berry phases. The latter has been verified in several 



theoretical and experimental applications in optics [|16H l8|l and elsewhere ||23|| . From these 
phase transformations U(l) of different components of the system we are able to obtain 
with proper combinations any desired U(2 m ) transformation. Since there exist experimental 
measurements of the Berry phase, it is plausible to expect the implementation of the U(2 m ) 
holonomic transformations. 

A further final advantage of the holonomic setup is that it is confined in the degenerate 
eigenspace produced by {|0), |1)}, describing one qubit. Entanglement of these states with 
the non-degenerate ones in the course of application of the logical gates does not occur due 
to the adiabaticity requirement. The initial control operators we use here, D, S, M and iV 
in general mix all the states of the Fock space, but at the end of the loop, only rotations 
between the degenerate eigenstates will be accounted for. 

The possibility to observe the proposed holonomies in the laboratory or even perform 
specific logical gates is a demanding task and an open question for the future. 
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